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ABSTRACT 

The Sunyaev-Zel'dovich effect is a powerful new tool for finding and studying clusters 
at high redshift, particularly in combination with their X-ray properties. In this paper 
we quantify the expected scaling relations between these properties using numerical 
simulations with various models for heating and cooling of the cluster gas. For a 
Non-radiative model, we find scaling relations in good agreement with self-similar 
predictions: Y oc T-^^^ and Y oc L^'^. Our main results focus on predictions from 
Cooling and Preheating simulations, shown by Muanwong et al. (2002) to provide a 
good match to the X-ray scaling relations at z = 0. For these runs we find slopes 
of approximately Y oc and Y oc Lx, steeper and fiattcr than the self-similar 
scalings respectively. We also study the redshift evolution of the scaling relations and 
find the slopes show no evidence of evolution out to redshifts well beyond one, while 
the normalizations of relations between the SZ signal and X-ray properties do show 
evolution relative to that expected from self-similarity, particularly at z < 1. 

Key words: hydrodynamics - methods: numerical - X-rays: galaxies: clusters - cosmic 
microwave background 



1 INTRODUCTION 

The prospect of surveying for galaxy clusters across a wide 
range of redshifts using the Sunyaev-Zel'dovich (SZ) effect 
(Sunyaev & Zel'dovich 1972, 1980) raises the hope of a much- 
improved understanding of cluster physics, and in particular 
of the evolution of the cluster gas. The SZ effect is a par- 
ticularly powerful tool when combined with observations of 
X-ray emission from the hot cluster gas, as the two tech- 
niques probe different properties of the cluster gas distribu- 
tion, and with the XMM-Newton and Chandra X-ray satel- 
lites in their main operations phases the number of clusters 
and groups with quality X-ray observations is rapidly in- 
creasing. It is therefore timely to extend theoretical models 
of clusters, where efforts have hitherto been concentrated 
upon their X-ray properties, to incorporate predictions for 
the SZ effect. 

Numerical simulations indicate that the distribution of 
non-baryonic dark matter in clusters, which dominates their 
mass, is approximately self-similar (e.g. Navarro, Frenk & 
White 1995, 1997). Observations indicate, however, that the 
baryonic gas component cannot share such a degree of self- 
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similarity. This is particularly evident from observations of 
the Lx-Tx relation, which show it to be steeper than pre- 
dicted by the self-similar model (e.g. Edge & Stewart 1991; 
Xue & Wu 2000). The reason for this discrepancy is that 
the gas is less centrally concentrated than the dark mat- 
ter, due to physical processes (in addition to gravity) that 
raised the entropy of the gas (e.g. Evrard & Henry 1991; 
Kaiser 1991; Bower 1997; Voit et al. 2002; Ponman, Sander- 
son & Finoguenov 2003). This entropy could come from di- 
rect heating of the hot gas from stars or AGN (e.g. Wu, 
Fabian & Nulsen 2000; Bower et al. 2001; Quihs, Bower & 
Balogh 2001), or from radiative cooling which transforms 
low-entropy material into stars allowing high-entropy mate- 
rial to flow in to replace it (e.g. Knight & Ponman 1997; 
Pearce et al. 2000; Bryan 2000). Understanding the relative 
importance of heating and cooling, and how these processes 
conspire to establish excess entropy in clusters, is one of the 
key problems in cluster gas physics. 

The implications of an excess of entropy on the global 
properties of the SZ effect (e.g. source counts, the mean 
Compton parameter and the angular power spectrum) has 
previously been investigated using both analytical tech- 
niques (e.g. Cavaliere & Menci 2001; Holder & Carlstrom 
2001; Zhang & Wu 2003) and numerical simulations (e.g. 
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Springel, White & Hernquist 2001; da Silva et al. 2001; 
White, Hernquist & Springel 2002). The mass dependence 
of the SZ effect in clusters has been addressed by Metzler 
(1998) and by White et al. (2002). More recently, McCarthy 
et al. (2003a) used a semi-analytic model to derive scaling re- 
lations between the central Compton parameter and mass, 
temperature, SZ flux density and X-ray luminosity. These 
scalings were subsequently used in a companion paper (Mc- 
Carthy et al. 2003b) to estimate the level of the entropy 
floor compatible with present-day SZ observations. 

The purpose of this paper is to investigate the rela- 
tionship between SZ and X-ray properties of clusters in A'^- 
body/hydrodynamical simulations that include models for 
both cooling and preheating. The simulations we use have 
already been demonstrated to give good overall agreement 
with observed X-ray scaling relations at redshift zero (Muan- 
wong et al. 2001; Thomas et al. 2002; Muanwong et al. 2002). 
Our approach is mainly focused on the derivation of theoret- 
ical scaling laws, and their evolution with redshift, between 
the integrated SZ flux density and other cluster properties, 
such as X-ray temperature and luminosity. We quantify de- 
viations from self-similar evolution on both the slope and 
normalization of the scalings and provide flts which can be 
used for comparison with other theoretical models and ob- 
servations. An observationally- motivated analysis of SZ/X- 
ray correlations, that can be more readily applied to 'blind' 
searches for high-redshift clusters, will be presented in a fu- 
ture paper. 

This paper is organized as follows. We define the ob- 
servational quantities used and the simple scaling laws pre- 
dicted by the self-similar model in Section 2. Details of our 
simulations and how our cluster catalogues were constructed 
are presented in Section 3. In Section 4, we investigate the 
correlation between thermal SZ integrated fiux and other 
3D cluster properties (mass, mass-weighted and emission- 
weighted temperature and X-ray luminosity) from the sim- 
ulation data at redshift zero. The evolution of these relations 
with redshift is studied in Section 5, before we draw conclu- 
sions in Section 6. 



2 THEORETICAL FRAMEWORK 

2.1 Definitions of physical quantities 

The SZ effect arises due to inverse Compton scattering of 
CMB photons off free electrons, with the largest SZ signal 
being produced by intracluster gas. The total SZ flux density 
produced by a cluster is the integral of the SZ sky brightness 
over its subtended solid angle 



S, = Io / [g{x)y-h{x)b] dfl, 



(1) 



where lo = 2(fcBTo)'^ /{hcf ~ 2.28 x 10* mjy arcmin"^, ^ 
for a mean CMB temperature of To ~ 2.725 K (Mather et 
al. 1999). The first term inside the brackets accounts for the 
thermal SZ effect, due to internal motion of the electrons, 
whereas the second term, the kinetic SZ effect, is the contri- 
bution due to bulk motion of the gas. We ignore the kinetic 



SZ effect in this paper as the thermal SZ effect dominates at 
all frequencies except v ~ 217GHz {x — hv/k^Ta ~ 3.83). 

The Comptonization parameter y contains information 
on the structure of the intracluster gas 



y = 



r„ n„ dL 



(2) 



where To and ric are the temperature and density of the 
electrons, ctt — 6.65 x 10~^^ cm'^ the Thomson cross-section, 
c the speed of light and the electron rest mass. It is 
therefore convenient to define the integrated contribution to 
5*1. from y 

kso-T 



Y 



y dO, - 
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d\ 



dA- 



TcTic dV, 



(3) 



where is the angular diameter distance and dA — df2 d\ 
is the sky projected area element of the cluster. Since Y 
depends on distance, we will usually work with the intrinsic 
thermal SZ signal, defined as F'"* = Y d\{z). 
The X-ray luminosity of a cluster is 



: J n,nuA,{T)dV. 



(4) 



The integral is over the volume of the cluster and Ac(r) 
denotes the cooling rate of the gas, assuming coUisional ion- 
ization equilibrium. At temperatures above a few keV, the 
emission is dominated by thermal bremsstrahlung, but be- 
low this line emission becomes important. 

2.2 Self-similar scaling relations 

In the absence of any non-gravitational heating and cooling 
processes, clusters scale self-similarly to a good approxima- 
tion (Kaiser 1986; Navarro et al. 1995). In self-similar mod- 
els, the temperature of the gas scales with the cluster mass 
as 



(5) 



assuming the density of the gas scales with the mean density 
p oc (l + z)^ and the system is in virial equilibrium. Applying 
this to Eq. @ we get 



ymt 
OC 



/gas r5/2 {l + z 
/gasM5/3(l-f , 



-3/2 



(6) 



where /gas is the gas mass fraction of the cluster. 

The X-ray luminosity of a self-similar cluster scales as 



3/2 



./glsT2(l + z) 

flsM^ri{i + z)y\ 



(7) 



1 mJy : 



10-2"^ erg s" 



assuming that the X-rays are due to bremsstrahlung emis- 
sion (Ac oc r^/2). 



3 METHOD 

3.1 Simulation details 

Results are presented from three simulations, described in 
detail elsewhere (Thomas et al. 2002; Muanwong et al. 2002); 
we summarize pertinent details here. 

We assume a ACDM cosmology, setting the density pa- 
rameter, rim = 0.35, cosmological constant, = 0.65, hub- 
ble parameter, h = 0.71, baryon density, Q,bh^ = 0.019, 
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CDM power spectrum shape parameter, F = 0.21, and 
normalization ug = 0.9. These values are similar to those 
favoured by current observations including WMAP. Initial 
conditions were generated using 160^ particles each of bary- 
onic and dark matter, perturbed from a regular grid of co- 
moving length, L = 100/i~^Mpc. These choices set the gas 
and dark matter particle masses equal to 2.6 x 10^?i~'^Mq 
and 2.1 x 10^°/i"^Mo respectively. 

The runs were started at redshift z = 49 and evolved 
to z = using a parallel version of the hydra N- 
body /hydrodynamics code (Couchman, Thomas & Pearce 
1995; Pearce & Couchman 1997), which uses a Smoothed 
Particle Hydrodynamical (SPH) implementation similar to 
that discussed in Thacker & Couchman (2000), and con- 
serves both energy and entropy to a satisfactory degree. The 
gravitational softening was fixed at 50/i~^kpc in comoving 
co-ordinates until z = 1, then at 25/i~^kpc in physical co- 
ordinates until z = 0. 

The simulations differed only in the adopted heat- 
ing/cooling model. In the first model, a Non-radiative sim- 
ulation, the gas was subject only to gravitational, adiabatic 
and viscous forces, and so could only change its entropy 
through shock-heating. This model fails to reproduce the 
X-ray scaling relations (Muanwong et al. 2001) but never- 
theless can be used to test the self-similar scalings presented 
in the previous section. 

The second simulation, a Cooling model, allowed gas 
to cool radiatively as described by Thomas & Couchman 
(1992). Tabulated values of Ac were generated from cool- 
ing tables in Sutherland & Dopita (1993), assuming a time- 
varying global gas metallicity, Z — 0.3{t/to)ZQ, where 
to ~ 12.8 Gyr is the current age of the Universe. Cooled 
material (identified as all gas particles with temperatures 
T < 1.2 X 10" K and overdensities 5 > 1000) was converted 
into coUisionless particles. Note that the choice of density 
threshold (considerably lower than the observed density of 
star formation) largely circumvents the problem of unphys- 
ical cooling of hot gas particles in contact with cold dense 
gas, as well as increasing run-time efficiency. The global frac- 
tion of cooled gas in the simulation box is ~15 per cent at 
z — 0. The removal of this gas from the hot phase sufficiently 
increases the entropy in the centres of clusters to reproduce 
the X-ray scaling relations, but the cooled fraction is 2 to 
3 times higher than observed (Cole et al. 2001). In this re- 
spect, this run provides a maximum variation of Y from a 
cluster due to effect of radiative cooling alone. 

Finally, the third simulation, a Preheating model, also 
allowed the gas to cool radiatively, but the gas was impul- 
sively heated by 1.5 kcV per particle at z — 4. This gener- 
ates the required minimum entropy in the centres of clusters 
but prevents further cooling to 10*K; the global fraction of 
cooled gas in the simulation box is only ~0.5 per cent at 
z = 0. 

Note that we do not attempt to perform a numerical 
resolution study for this work. Apart from being compu- 
tationally demanding, it would be misleading to study the 
effects of resolution on our current models. In particular, 
increasing the resolution in the Cooling model would result 
in an even higher cooled fraction, and so attaining conver- 
gence in this model is undesirable. Instead we interpret our 
models as a relatively crude first attempt at capturing the 
effects of increasing the entropy of the gas by the required 



amount in two completely different ways (i.e. from cooling 
and heating). A proper convergence analysis is warranted 
when feedback is included (e.g. Springel & Hernquist 2003) 
and will be investigated in future work. 



3.2 Cluster selection and SZ/X-ray estimators 

We identified clusters in our simulations according to the 
method described in Thomas et al (1998). Clusters were 
selected by first creating a minimal-spanning tree of all dark 
matter particles whose density exceeds S = 178^""'^'^ (z) 
times the mean dark matter density (the approximate value 
for a virialized sphere, as predicted by the spherical top- 
hat model; Eke, Navarro & Frenk 1998). The tree was then 
pruned into clumps using a maximum linking length equal to 
0.55^^^'^ times the mean intcrparticlc separation. A sphere 
was then grown around the densest particle in each clump 
until the enclosed mean density exceeded a value A in units 
of the comoving critical density. For this paper, we use A = 
200, larger than the virial value for our cosmology (A ~ 
111) but commonly used by other authors. Master cluster 
catalogues^ were produced containing only objects with at 
least 500 particles each of gas and dark matter, equivalent 
to a mass limit. Mum ~ 1-18 x lO"ft"^M0. At z = 0, the 
catalogues contain 428, 457 and 405 clusters for the Non- 
radiative, Cooling and Preheating simulations respectively. 

For each cluster, we then estimated various observable 
quantities: the intrinsic SZ signal, F'"*, the X-ray emission- 
weighted gas temperature, Tx and the X-ray luminosity, Lx. 
We identified the hot intracluster gas as all gas particles 
within i?2oo with T > 10^ K (in practice there are very few 
particles with 10'*K< T < lO^K as their cooling times are 
very short) and assume full ionization, such that the number 
of electrons per baryon is ?7 = 0.88 for a hydrogen mass 
fraction X = 0.76. We have 



int 



l^JL f MT 

o J gas -tmwj 



(8) 



meC^ rup 

where = ^ rniTij ^ nii is the mass- weighted tem- 

perature of the gas, /gas = J2 'mi/M the gas fraction, rrn 
the mass of hot gas particle i and M the total mass of the 
cluster. 

The X-ray emission-weighted temperature was esti- 
mated as 



Tx = 



"EimiPiAsoftiTi, Z)Ti 
EimiPiAsofl(Ti,^) ' 



(9) 



where p, and Tj are the density and temperature of the hot 
gas particles. Since most observed temperatures use instru- 
ments sensitive to soft X-rays, we use a soft-band cooling 
function, Agott, from Raymond & Smith (1977) for an en- 
ergy range 0.3-1.5 keV. 

The bolometric X-ray luminosity of each cluster, cor- 
rected from its soft-band emission, was estimated as 



Aboi(2x 



E 



miPiAsoit(Ti, Z) 



(10) 



Asoft(Tx) ^ {p-VTiYlY 

where prriH = 10~^'*g is the mean molecular mass of the gas 

^ The catalogues are available to download at the website 
Virgo . Sussex .ac.uk 



© 0000 RAS, MNRAS 000, 000-000 



4 A. C. da Silva et al. 



and Aboi is the bolometric cooling function (Sutherland & 
Dopita 1993). 

Finally, to avoid introducing artificial trends in the de- 
rived scalings with temperature and X-ray luminosity (Sec- 
tions 4 & 5), we imposed lower limits on Tx and Lx (by 
inspecting the Tx-M and Ly^-M relations at all redshifts) 
to create (separate) catalogues that were complete in the 
two quantities. At each redshift, these limits approximately 
correspond to the maximum temperature and luminosity of 
all clusters at our mass limit, A/jim. For example, when in- 
vestigating temperature scalings at z = 0, all clusters with 
Tx < 0.35, 0.65, 0.60 keV (for the Non-radiative, Cooling and 
Preheating simulations respectively) were discarded, reduc- 
ing the number of clusters in our catalogues by around 30 per 
cent in the Cooling and Preheating simulations and 46 per 
cent in the Non-radiative run. Similarly, for scalings with X- 
ray luminosity at 2 = 0, our original cluster catalogues were 
trimmed by selecting only clusters with Lx > 2.9 x 10*^, 
9.0 X 10""^, 4.5 X 10""^ erg s~^, for the Non-radiative, Cooling 
and Preheating simulations, reducing the number of clusters 
in the original catalogues by about 52, 46, and 29 per cent 
respectively. 



4 SCALING RELATIONS AT REDSHIFT ZERO 

4.1 The Y-M and Y-Tmw relations 

We begin by correlating F'"* with mass and mass-weighted 
temperature, all measured within R2oo- An early study of 
the relation between mass and SZ absorption in clusters was 
made in unpublished work by Metzler (1998), who verified 
with non-radiative simulations the related self-similar scal- 
ing, y oc M. More recently. White et al. (2002) also recov- 
ered the self-similar scaling, Y cx M'^^^ , from simulations 
with the same box size but a slightly higher (factor of 2.5) 
mass resolution than that used here, concluding that the 
combined effect of cooling and feedback (galactic winds) on 
the SZ properties of their clusters was small. 

Fig.Qillustrates our F200 ~ Af2oo and Y — Tmw relations 
for clusters in the Non-radiative (crosses) and Cooling (cir- 
cles) simulations. The solid lines represent power-law best 
fits to the data, while the dashed line is the best-fit to clus- 
ters in the Preheating simulation, omitted for the purpose 
of clarity. The figure shows that Y^ is tightly correlated 
with mass and mass- weighted temperature in all runs. As we 
shall see, however, these correlations show significantly less 
scatter than those relating Y with X-ray emission-weighted 
properties, because the latter is sensitive to substructure in 
the dense core gas. The best power-law fits to the y200~-^200 
relation in each simulation are 

• Non-radiative run: 
Y^ = 1.22 X 10-« ( V''' {h-'Mpcf, (11) 



1014 /1-IM0 



Cooling run: 



Y^o = 0.97 X 10-« ( ^ {h-' Mpc) ^ (12) 



^0 



Preheating run: 



YiS'o = 1.12 X 10"'' 
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Figure 1. Scaling relations between V2130 ' 



M200 (top panel) and 



¥"200 ^ Tinw (bottom panel) for the Non-radiative (crosses) and 
Cooling (circles) simulations. Solid lines are power-law fits to the 
data from these runs, while the dashed line is the corresponding 
fit for clusters in the Preheating simulation. 



As expected, the slope of the fit to the Non-radiative 
clusters is close to 5/3, the self-similar value given in Eq. Q. 
The slope is steeper for the Cooling and Preheating simula- 
tions however, primarily due to lower Y values for low-mass 
systems. This difference is due to the effects of cooling and 
heating on both the gas fraction and temperature of the 
clusters. Cooling removes low-entropy (core) gas from the 
halo, causing higher-entropy material to flow in to replace 
it. This results in a lower gas fraction although the remaining 
gas is hotter (Pearce et al. 2000; Bryan 2000; Muanwong et 
al. 2001). For low-mass clusters (M200 < lO^^/i'^M©), the 
net effect of cooling on Y is dominated by the decrease in hot 
gas fraction, shown in Fig. \^ around a factor of two lower 
than the Non-radiative average at M200 ~ 10^^ h~^ Mq. For 
the highest mass clusters this difference decreases to around 
30 per cent; combined with the increase in temperature of 
the gas, the clusters have very similar Y values in the two 
simulations. 

In the Preheating case, the energy injection at high red- 
shift heats the gas and expels some of it from the cluster. 
Hence its effect on the cluster properties is similar to the 
cooling model but the fate of the gas is different (Muan- 
wong et al. 2002). Fig. |21 shows that the gas fraction in Pre- 
heating clusters increases more rapidly with mass than for 
objects in the Cooling simulation (please refer to Section 
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Figure 2. Fraction of hot gas within -R200 a function of mass 
for clusters in the Non-radiative (crosses), Cooling (circles) and 
Preheating (diamonds) simulations. 



3.1 and Fig. 3 of Muanwong et al. 2002 for a detailed dis- 
cussion of the hot, cold and total baryonic mass fractions 
in clusters for the same set of simulations). In low-mass 
systems the heating energy was sufficient to significantly 
decrease /gas (and hence Y), but left gas fractions almost 
unchanged in the highest-mass clusters, where the heating 
resulted only in a modest increase in temperature. Note, 
therefore, that the best-fit relation is only a good descrip- 
tion of the low-mass clusters in the Preheating simulation; 
systems more massive than ~ lO^'*/i~^M0 are adequately 
described by the Non-radiative relation. This also explains 
why White et al. (2002) find good agreement with the self- 
similar Yint — M scaling, as they only investigated high-mass 
clusters {M200 ^ 5 x 10 H'^Mq), where a high fraction of 
the wind material would have been retained. 

The efi'ects of cooling and preheating on the cluster 
properties are more readily apparent in the relation be- 
tween Y and mass-weighted temperature. For a cluster of 
fixed gas and total mass, Y oc Tmw, and so heating the gas 
causes a shift to the right of the steeper self-similar rela- 
tion (F oc T^''^). Additionally, the decrease in gas fraction 
lowers Y, separating the two relations further. The best-fit 
^200 ~ Tmv relations are 



Non-radiative run: 



Knn = 5.59 x 10" 



• Cooling run: 
YiSo = 2.51 X 10-' 

• Preheating run: 
K,™ = 2.23 X 10-" 



IkeV 



IkeV 



IkeV 



(/i-'Mpc)^ 



(14) 



(15) 



(16) 



4.2 The Y -Tx relation 

Of more practical interest is to investigate how Y corre- 
lates with X-ray observables, particularly luminosity and 
emission-weighted temperature. 

Fig. 13 illustrates the relationship between V200 aud Tx 
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10 



Figure 3. Fjqq versus Tx for clusters in the Cooling (top panel) 
and Preheating (bottom panel) simulations, including (circles) 
and excluding (crosses) soft-band emission from within the cool- 
ing radius. Also shown for comparison is the scaling with the 
mass-weighted temperature found in the Non-radiative simula- 
tion (triangles). 



for the Cooling (top panel) and Preheating (bottom panel) 
simulations. Circles represent clusters for which all hot gas 
particles within 7?200 are considered, whereas crosses are for 
the same clusters but excluding particles within the cooling 
radius in the computation of the temperature. We define 
the cooling radius as the radius within which the gas has a 
mean cooling time of 6 Gyr. (This procedure is an attempt 
to provide a 'cooling flow' correction to the estimated values 
of Tx.) 

Comparing these results with Fig. Q we see that the 
scatter is larger for the correlation with Tx than with Tmw- 
This is because Tx is weighted by density and temperature 
and is therefore more sensitive to substructure, particularly 
from the centres of clusters. This is also evident from the fact 
that the scatter is smaller when the gas within the cooling 
radius is excluded (note that this also increases Tx in high- 
mass systems, which slightly flattens the Y — Tx relation). 

Table lists best-fit coefficients, ut and At, to the 
power-law, y200 = (At/10") (fcBTx/lkeV)"T for the Cooling 
and Preheating runs. (We also present results from the Non- 
Radiative run, using a bolometric cooling function, and this 
agrees well with the self-similar scaling Y oc T^^'^.) For a 
given run, the normalization of the y200 ~ 2x relation is 
lower and slope steeper than for the V200 ^Tmw relation. The 
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Table 1. Power-law fits to the simulated cluster Y — Tx and 
Y — Lx scalings at ^ = 0. Here and ol are the slopes of 
the VjQQ — Tx and V^'qq — Lx relations, whereas and 
Al/IO'^ are the values of Yj'oc 
/cbTx = IkeV and Lx = W*^h 



^ erg s ^ , respectively. 







^int 


-Tx 






Run 




Ox 


At 


«L 




Non-radiative 




2.47 


10.1 


1.26 


0.37 


(bolometric) 












Cooling 


uncorrected 


3.35 


1.64 


0.94 


4.72 


(soft-band) 


cooling-flow 












corrected 


3.06 


1.41 


1.09 


6.51 


Preheating 


uncorrected 


3.26 


1.93 


0.98 


6.65 


(soft-band) 


cooling-flow 












corrected 


3.18 


1.70 


1.04 


8.09 



heavier weighting of the denser gas exacerbates the eflects 
of cooling/heating on the Y — T relation discussed in the 
previous section. 



4.3 The y - Lx relation 

We end this section by reporting on the z = correlation be- 
tween y'"* and bolometric X-ray luminosity, Lx, as inferred 
from the soft-band emission. Fig. m shows the yjoo - Lx 
relation from the Cooling (top panel) and Preheating (bot- 
tom panel) simulations. Again, we use circles to represent 
values including all hot gas particles, whereas crosses rep- 
resent results when excluding gas from within the cool- 
ing radius. The best-fit parameters and j4l, defined by 
YiSo = (^l/10") (Lx/L43)°^ where L43 = lO*^/*-" ergs-\ 
are listed in Table Q 

The Cooling and Preheating relations are significantly 
shallower than the predicted self-similar scaling, Y oc L^''*, 
which is well reproduced by the Non-radiative simulation 
(where ql ~ 1.26). Again, this is due to the increase in en- 
tropy of the gas, particularly in low-mass clusters, which 
decreases their luminosity. Omitting gas from within the 
cooling radius slightly increases the slope of the Y'"' — Lx 
relation and also reduces the scatter. 



5 EVOLUTION OF SCALING RELATIONS 

We now investigate the evolution of the scaling relations 
with redshift. In each of the following subsections, we fac- 
tor out the dependence expected from self-similar evolution. 
Note that the self-similar evolution has no dependence on 
fio because we define the overdensity with respect to the 
comoving critical density. 



10" 



10 



>- 10 



10 



cooling 



+ Of o 



totol soft-bond emission o 
exclude emission within R„qi + 



10^ 



10-^ 10" 

r,40 , -2 



10- 



Lx / (10"° h-' erg s" 



10 



10 



>- 10 



10 



preheat. 



+ 40 

©4-0 




totol soft-bond emission o 
exclude emission within Rcoq, + 



10-^ 10"* 
/ (10"° h-^ erg s"') 



10- 



Figure 4. V200 versus Lx for clusters at 2 = in the Cooling 
(top panel) and Preheating (bottom panel) simulations. Circles 
(crosses) denote objects when gas within the cooling radius is 
included (excluded) in the calculations. 
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5.1 Evolution of the Y — M relation 

Fig. shows the Yioo — M200 relation for clusters in the 
Cooling simulations, for a range of redshifts between z ~ 
and z ^ \. Here we use "^200 = ^200^X^5 demonstrating the 
dominant effect of dA on Y at these redshifts. Note that 



Figure 5. Y200— relation at various redshifts for the Cooling 
simulation. The solid line is the power-law best fit. 
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Figure 6. Redshift dependence of the slope (om) and normaliza- 
tion {/3m) of the cluster 5^200 "-'^200 relation in the Non-radiative 
(crosses), Cooling (circles) and Preheating (diamonds) simula- 
tions. Symbols represent best-fit parameters and the dotted lines 
their 1-a uncertainties (see text). The dashed line represents the 
slope predicted by the self-similar model (ctM = 5/3). 



Table 2. Deviations from self-similarity in the mean slope and 
normalization of the SZ scaling relations. The quantities (aM)i 
(ax) and (ol) are the mean values of the best- fit slopes of the 



200 



M200, y. 



200 



; and 5^200 ~ relations respectively, 
averaged over redshifts 2 = to z = 2. The quantities tm, 7t 
and 7l are power-law best fits to the normalization of the scalings, 
10'3(^) (X (f + 2)^; the first two are a good fit for the complete 
redshift range, but for luminosity the fit is restricted to < 2 < 1. 



Run 


Mean slope deviation from self-similarity 
(om) - 5/3 (ot) - 5/2 («l) - 5/4 


Cooling 
Preheating 


0.f6 0.44 -0.16 
0.28 0.52 -0.22 


Run 


Normalization deviation from self-similarity 
7m 7t 7l 


Cooling 
Preheating 


0.14 0.73 1.30 
0.02 0.33 1.62 



more apparent if the SZ flux density is evaluated within a 
smaller radius, where the average gas density is higher. As 
we shall see, such effects are apparent when studying evolu- 
tion of the X-ray properties of the clusters. 



the number of clusters satisfying our mass selection crite- 
ria, A/200 1-18 X W^'^h~^ Mq, decreases with increasing 
redshift. 

A noticeable feature of Fig. I^is the apparent indepen- 
dence of the slope of the 1200 — M200 relation with redshift. 
To quantify the effects of cooling and heating on the evolu- 
tion of the slope, we flt a power-law 



-'^200 
1 + z 



10 



/3m (z) 



(M2oo/lO''*/i"'Af0) 



(17) 



to the distribution of clusters at each redshift, and plot the 
resulting best-fit values of au against z in the upper panel 
of Fig. |n] (We do not plot results for redshifts beyond 2 as 
the number of clusters below our mass threshold becomes 
too small for a reliable fit to be produced.) Crosses, circles 
and diamonds represent best-fit values, and dotted lines de- 
note 1-0" uncertainties (estimated using the method given 
by Press et al. 1992) from the Non-radiative, Cooling and 
Preheating simulations respectively. Indeed, the results are 
consistent with an unevolving slope for all three models out 
to 2 = 2. 

The lower panel of Fig.|n]shows the evolution of the nor- 
malization, obtained by fitting the same power-law function 
as before, but with a fixed slope taken to be equal to the 
mean value of ayiiz) within the displayed range of z (val- 
ues are given in Table From this figure, we see that the 
evolution of the ^200 — Af200 normalization is remarkably self- 
similar for all runs. Fitting a power-law, 10''"'^' cx (1-1- 2;)™, 
to the derived normalizations (within the displayed range of 
z), we find that the largest deviation from self-similar evolu- 
tion of the normalization is only 7m ~ 0.14, for the Cooling 
simulation (Table I^J. 

It is therefore clear that the evolution of the lint — 
M200 relation is not a sensitive probe of underlying non- 
gravitational physics, as it is a measure of the large- 
scale evolution of clusters. As discussed in McCarthy et 
al. (2003a), the effects of cooling and preheating are much 



5.2 Evolution of the Y — Tx relation 

Fig. 13 shows the redshift dependence of the slope and nor- 
malization of the Y — Tx relation for the Cooling and Pre- 
heating simulations, using cooling-flow corrected tempera- 
tures. As in the previous section, we flrst flt the cluster dis- 
tribution at each redshift with a power-law 



(l + z)-3/2 



= IO^t(^) (Tx/keV)"^<"> , 



(18) 



to determine the evolution of the slope. While this shows 
signiflcant statistical oscillations (reflecting the higher sen- 
sitivity of the emission-weighted temperature to the thermal 
history of the gas), there is no evidence of evolution with red- 
shift beyond that predicted by self-similarity. The slope is 
then fixed at the mean value, (qt), before determining the 
redshift dependence of the normalization, Pt{z). Unlike the 
y^'-Af relation, there is a deflnite evolutionary trend that is 
particularly signiflcant in the Cooling case. This is an indica- 
tion that evolution is predominantly in the X-ray properties 
rather then in the SZ signal, which is much less sensitive to 
the details of the gas distribution in the central regions of 
clusters. Fitting the normalizations with a power-law of the 
form 10''t(^) oc (1 z^^ , we find 7t = 0.73 and 7t = 0.33 
for the Cooling and Preheating simulations, respectively. In 
both cases the normalization of the Tx-M relation increases 
with time compared to the expected self-similar evolution. 



5.3 Evolution of the Y — Lx relation 

We end by studying the evolution of the Y — Lx relation, 
shown in Fig. |H| Our power-law fltting function was 



V" 



(1 + ^) 



-27/8 



10'3^(^' (Lx/L43)"^^^' , 



(19) 



using cooling-flow corrected values for Lx- Again, the slopes 
are constant with redshift to good approximation, and there 
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Figure 7. Rcdshift dependence of the slope (ax) and normal- 
ization (/3t) of the cluster Y^qq — Tx relation in the Cooling 
(circles) and Preheating (diamonds) simulations. Symbols repre- 
sent best-fit values and dotted lines their f-cr uncertainties. The 
dashed line represent the slope predicted by the self-similar model 
(ax = 5/2). 



Figure 8. Redshift dependence of the slope (cvl) and normal- 
ization (/3l) of the cluster Y2qq-Lx relation in the Cooling (cir- 
cles) and Preheating (diamonds) simulations. Symbols represent 
best-fit parameters and dotted lines their l-cr uncertainties. The 
dashed line represent the slope predicted by the self-similar model 
(aL = 5/4). 



is additional evolution in the normalization. In this case a 
single power-law is not so good a fit over the complete red- 
shift range, and we restrict the fit to < 2; < 1 which is 
where observations are likely to be made. Fitting a power- 
law lO^Lt^) (x{l + zY'-^ yields 7l ~ 1.3 & 1.6 in the Coolmg 
and Preheating cases respectively (Table |5J in this redshift 
range, with the relation flattening at redshifts beyond one. 



6 CONCLUSIONS 

In this paper, we investigated the relationship between SZ 
and X-ray properties of clusters drawn from three simula- 
tions: a Non-radiative simulation, where the gas could only 
heat adiabatically and through shocks; a Cooiinj simulation, 
where the gas could also cool radiatively and a Preheating 
simulation, where we additionally heated the gas by 1.5 keV 
per particle at 2 = 4. We focussed on the relations between 
the thermal SZ signal Yjnt and cluster mass, temperature 
and X-ray luminosity, both at z = and the evolution with 
redshift. Our main conclusions are as follows: 

• Yint is tightly correlated with both mass and mass- 
weighted temperature. Scaling relations of Vint with X-ray 
quantities show significantly more scatter however, due to 
the sensitivity of the latter (which are emission-weighted) 
to the gas distribution in the cores of clusters. Excluding 
gas from within the cooling radius of each cluster reduces 
the scatter. 

• The Non-radiative simulation scaling relations closely 
match those predicted by the self-similar model, with 
Yint cc T^'^ and Vint oc L^'^ . 

• Cooling decreases Yint for a given cluster mass, due to 
the reduction in the hot gas fraction (which decreases Yint) 
being more effective than the increase in temperature of the 
remaining gas (which increases lint). Preheating the gas has 
a similar effect except in the most massive clusters, where 
the heating energy was not sufficient to expel much gas from 



the halo. At redshift zero, the resulting Yint-Tx scaling rela- 
tion is steeper than the self-similar relation, Yint oc T^ °~^'^, 
and the Yint-ix relation flatter, Kint c>c L^®~^'^. 

• There is no evidence of evolution in the slopes of any of 
the scaling relations. The normalization of the 1200 — M200 
relation is also consistent with self-similar evolution. How- 
ever, the relations between Fint and the cluster X-ray prop- 
erties do show significant evolution in normalization relative 
to that expected from self-similarity, indicating that evolu- 
tion is predominantly in the X-ray properties rather than 
the SZ properties. The effects of excess entropy reduces the 
negative evolution of the Yint-7x and Yint-i'X with redshift, 
particularly at 2; < 1. 
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